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The filter response of nonuniform, almost-periodic structures, such as 
corrugated optical waveguides, is investigated theoretically. The filter 
process, leading to reflection of a band of frequencies near the Bragg fre- 
quency, is treated as a contradirectional coupled-wave interaction and 
shown to obey a Riccati differential equation. The nonuniformity of the 
structure is represented by a tapering in the coupling strength (e.g., the 
depth of the corrugation) and by a chirp in the period of the structure. 
For small reflectivities, the filter response is a Fourier transform of the 
taper function. For large reflectivities, the Riccati equation was evaluated 
numerically and plots are given for the response of filters with linear and 
quadratic tapers and ivith linear and quadratic chirps. 

I. INTRODUCTION 

Recent papers 1 - 2 have reported the fabrication of wavelength-selec- 
tive filters in the form of corrugated, thin-film, optical waveguides and 
have described the evaluation of the filter characteristics by means of 
tunable dye lasers. These papers include a numerical comparison of 
experimental and theoretical results on the nature of the filter response 
caused by a gradual tapering in the corrugation or grating strength 
(which we call "tapers"), or by a gradual variation of the effective 
grating period (called "chirping"). The present paper describes in some 
detail the theory by which these numerical results were obtained; it 
offers a discussion of the general characteristics of nonuniform grating 
filters, such as their scaling properties and the symmetry of the filter 
response ; and it provides a more complete collection of plots represent- 
ing numerical results for linear and quadratic tapers in the coupling 
strengths as well as linear and quadratic chirps in the effective grating 
period. Periodic waveguides serve as band-rejection filters with a 
center frequency determined by the Bragg condition 

K = 20(c), (1) 
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where K = 2ir/A is the grating constant or spatial frequency, A the 
grating period, and /3(w) the propagation constant of the waveguide 
mode of interest. The filter mechanism results from the backward 
scattering of the light by the periodic structure which we describe as a 
contradirectional coupled-wave interaction. For its analysis, we use 
the coupled-wave formalism 3 ' 4 which allows adequately for the deple- 
tion of the incident light. 

Problems similar to ours have been encountered in other coupled- 
wave processes. Codirectional coupling in nonuniform structures has 
been investigated in connection with microwave directional couplers, 5-7 
holography, 8 and tapered optical directional couplers. 9 Our present 
problem of contradirectional coupling in nonuniform periodic structures 
has been considered by Uchida 10 who analyzed structures with ex- 
ponential tapers and by Hill and coworkers 11-13 who use a numerical 
method based on the iteration of a pair of coupled-mode integral 
equations. Another possible method is to approximate the nonuniform 
structure by a set of short uniform elements, each represented by a 
known matrix, and use matrix multiplication to calculate the proper- 
ties of the compound overall structure. Kermisch 8 has applied such a 
technique to codirectional coupling in hologram gratings. As described 
below, we have used yet another method where we reduce the pair of 
coupled-wave equations to a single Riccati differential equation, which 
can then be solved by tested numerical techniques such as the Runge- 
Kutta method. 

II. RICCATI EQUATION FOR THE REFLECTION COEFFICIENT 

Following coupled-wave theory, we assume a periodic, single-mode 
waveguide with an electromagnetic field which can be represented by 
two contradirectional coupled waves in the form 

E(z) = R(z)e-*<>* + S(z)e#°°, (2) 

where R and S are the complex amplitudes of the forward- and back- 
ward-running mode. These amplitudes are linked by the standard 
coupled- wave equations 3,4 

R' + j8R jkS(T>+ (3) 

and 

S' - J8S = JKRe», (4) 

where we have allowed for a phase-shift of the periodicity relative to 
the origin (z = 0). The prime indicates differentiation with respect to 
2. The coupling coefficient k is related to the amplitude of the wave- 
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guide perturbation (e.g., the height of the film corrugation), and 
formulas for specific guides with specific perturbations are given in the 
literature. 3,4,14 The measure 5(u>) indicates the frequency deviation 
from the Bragg condition and is defined as 4 

5 = j8 - K/2 = - /3o « Aco/« ff , (5) 

where /3o = K/2 is the propagation constant at the center (Bragg) 
frequency, Aw is the radian frequency deviation from that frequency, 
and v a is the group velocity of the guide. 

We consider, now, structures in which both the coupling coefficient 
k(z) and the grating phase <f>(z) are slowly varying functions of z, 
indicating the nonuniformity in the grating parameters. We assume 
that the structure has a length L and extends from z = —L/2 to 
2 = L/2. The boundary conditions for our scattering problem are then 

R(-L/2) = 1, S(L/2) = 0. (6) 

The key to the reduction of the coupled-wave equations to a single 
differential equation is the definition of a local reflection coefficient p(z), 

*-!*-*■ (7) 

The z-derivative of this is 

»' = (l-f -*'!)-• » 

Combining the above expressions with (3) and (4), we obtain a Riccati 
differential equation for p which is of the form 

p'=j(25-<*>')p+;k(1+p 2 ). (9) 

The boundary condition for this equation follows from (6) as 

p(L/2) = 0. (10) 

Our quantity of interest is the reflection coefficient p( — L/2) of the 
entire structure, or the corresponding reflectivity pp*. What we call 
filter response is the dependence of this reflectivity on the frequency 
(or wavelength) of the light. Results for the response curves of specific 
taper functions are provided in later sections, but first we study some 
of their general properties. 

III. FOURIER-TRANSFORM RESPONSE AT LOW REFLECTIVITIES 

When the reflectivities are low, we expect a simple relationship 
between the taper function /<(z)e ; * (z) and the response function p(8). 
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This relationship follows from the assumption of an undepleted incident 
wave or the use of a first Born approximation. 1516 It is easily derived 
from the Riccati equation (9) by substituting a new variable a de- 
fined by 

p = ae m*-*\ (11) 

which obeys the differential equation 

a' = jdQr* 9 *-** + •V» to -*>J (12) 

When reflectivities are low, the term proportional to a 2 can be neglected, 
and we can integrate this equation with the result 

tr(-L/2) = -j I ds-/cOs)-e-" M -*>. (13) 

J-L/2 

This shows that, apart from phase factors, the response function p(5) 
is the Fourier transform of the taper function k(z) -e^ (,) . 

IV. SYMMETRY OF FILTER RESPONSE 

Let us now consider the conditions for which the filter response is 
symmetric, i.e., for which 

PP*(3) = PP *(-B). (14) 

To simplify our discussion, we assume that the filter structure can be 
described by a real and positive function n(z) and an arbitrary phase 
function 4>{z). Then the Fourier transform relation (13) predicts that 
the filter response is symmetric (a) if <p (z) = 0, or (b) if both k and <f> 
are symmetric, i.e., if k(z) = k( — z) and <t>{z) = <p(—z). Relation (13) 
is valid for small reflectivities only, but we get the same answer from 
the Riccati equation (9), which is valid for all reflectivities. To dem- 
onstrate this, we take the complex conjugate of (9) and write the 
result in the form 

(-V)' = /(-ai + *')(-p*) +*p + (-p*) 2 ]. (is) 

If we replace 6 with — 5 and <p' with — <t>', this differential equation for 
(— p*) is the same as the differential equation (9) for p, and (as the 
boundary conditions are also the same) it predicts the same reflectivity 
pp*. While replacing 5 with —8 simply means that we are looking 
at the other side of the center frequency, replacing <p' with — <£' means 
that we have replaced the original filter structure with another one. 
For a structure with symmetrical k(z) and <p{z), we have <p' = —tp'( — z), 
and the above replacement means that we have physically reversed the 
filter. Since the structure is lossless and reciprocal, the reflectivity is 
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Fig. 1 — Linear taper functions k(z) for the taper constants T = 0, 0.5, 1, 1.5, and 2. 

the same at the input and output port of the filter, and condition (b) is 
proven. Condition (a) for </>' = is easily proven by inspection of (15). 

V. SCALING OF FILTER RESPONSE 

In the general case, the Riccati equation (9) has to be evaluated 
numerically. To make the numerical results more broadly applicable, 
it is convenient to introduce the normalized quantities z/L, 8L, kL, and 
<f>'L. It follows from (9) that two filters of different length (labeled 1 
and 2) have the same reflectivity at the frequencies diLi = 82L2 if their 
taper functions obey the scaling rule 



Ki {i) Li = Ki (r 2 ) L2 



(16) 
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VI. TAPERED FILTERS 

Tapered filters are structures where the coupling coefficient k(z) 
varies along the length of the device. For simplicity, we consider here 
tapered filters with no chirp (0 = 0). At the center frequency (5 = 0), 
we can write the Riccati equation (9) of such a filter in the form 

dp 



1+p2 =>( 2 )^, 
which is easily integrated to yield for the reflection coefficient 

rL/2 



/LI2 
dzx(z). 
-L/2 



(17) 



(18) 



This is similar to a result found by Kermisch 8 for codirectional inter- 
actions in holograms. 



(a) 
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Fig. 2 — Frequency response of filters with a linear taper for the taper constants 
T = 0, 0.5, 1, 1.5, and 2 and values of (a) k L = tt/4, (b) kqL = ir/2, (c) kqL = 3jt/4, 
and (d) kqL = x. 
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Fig. 2 (continued). 
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6.1 Linear tapers 

In a linearly-tapered filter structure, the coupling coefficient k(z) 
varies as 

k = ko(1 + Tz/L), (19) 

where the constant T indicates the degree of the taper. Figure 1 shows 
a sketch of this taper function for five values of T including the cases 
of no tapering {T = 0) and full tapering (T = 2). From (18) and (19), 
we calculate the reflectivity of linearly-tapered filters at the center 
frequency (5 = 0) as 

PP * = tanh 2 Oc„L). (20) 

According to rule (a) of the previous section, the filter response of 
tapered filters is symmetric. Figures 2a through 2b show the filter 
response of linearly tapered filters for positive frequency deviations 
and five values of the taper constant T. The curves are the result of a 
numerical evaluation of (9) using the fourth-order Runge-Kutta 
method. Response curves are shown for values of k L = n-/4, ir/2, 
37r/4, and ir, and we notice a washing-out of the zeros in the response 
curve with increasing amounts of tapering. 

6.2 Quadratic tapers 

A filter with a quadratic taper is characterized by a coupling coeffi- 
cient k(z) that varies as 

K = Ko(l-^+Tz>/L*y (21) 

The particular form of this expression has been chosen to make the 
reflectivity at center frequency (5 = 0) independent of the constant 
T, and equal to 

pp* = tanh 2 (K L) (22) 

as calculated from (18). Figure 3 shows quadratic taper functions for 
five values of T. For positive T the curves are concave upward, and for 
negative valued T they are convex upward. Figures 4a through 4d 
show the numerically evaluated response curves for the same values 
of T and four typical values of k q L. We note that the concave upward 
tapers produce very high side-lobe levels in their response characteris- 
tics, while the side-lobe levels can be very low for the convex upward 
tapers. The side-lobe levels of the tapers with T = — 6 (emphasized 
by the thicker lines) are the lowest with about 2 percent in reflectivity. 

VII. CHIRPED FILTERS 

A variation of the grating period along the length of the filter is 
called a "chirp." Chirped filter characteristics may also be due to a 
variation of a waveguide parameter such as the refractive index or 
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Fig. 3 — Quadratic taper functions k(z) for the taper constants T = —6, —3, 0, 6, 
and 12. The curves intersect at the points z = ±L/2-v3 and k = k . 
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the guide height or width. We have chosen the phase function <j>(z) to 
represent all variations of this kind. We are, therefore, dealing with a 
perturbation proportional to cos (Kz + <f>) with constant spatial fre- 
quency K and a variable phase shift 0(z), which can also be viewed 
as a perturbation with a variable, or chirped, spatial frequency 
K + AK(z), where 

AK(z) = *'(«), (23) 



or as a perturbation of variable period A + AA(z), where 

A A (z) _ __ A-K(z) 
A K 



(24) 
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Fig. 4 — Frequency response of niters with a quadratic taper for the taper constants 
T = — 6, -3, 0, 6, and 12, and values of (a) k L = tt/4, (b) k L = tt/2, (c) k L 
= 3x/4, and (d) «oL = x. 
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Fig. 4 (continued). 
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A chirped grating period implies that the (local) Bragg frequency 
changes along z, which results in a general broadening of the filter 
response. 

7.1 Linear chirp 

A periodic structure with a linear chirp is described by a linearly 
varying spatial frequency which we write in the form 

AK (z) = 4>'{z) = 2Fz/D, (25) 

where the constant F is a measure for the degree of the chirp. The 



-F/2tt=0.5 



kL=tt/4 



(a) 
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NORMALIZED FREQUENCY 6L 



Fig. 5 — Frequency response of filters with a linear chirp for the chirp constants 
F/2* = 0, 0.5, 1, 1.5, 2, and 3, and for *L-values of (a) t/4, (b) w/2, (c) 3tt/4, and 
(d)m. 



120 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1976 



«L=3tt/4 




1 2 3 4 5 6 7 

NORMALIZED FREQUENCY 6L 



Fig. 5 (continued). 
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corresponding phase function is 

<f>(z) = Fzyu. 
For the grating period we have 

AA(g) = _ AK(z) _ _ F A 2z 
A K 2ir' LL ' 



(26) 



(27) 



For the difference AA of the grating periods at the center (z = 0) and 
end (z = L/2) of the device, we obtain 



AA^ A(0) - A (L/2) = F_ A 
A A 2ir'L' 



(28) 




-4 4 

NORMALIZED FREQUENCY 6L 



Fig. 6 — Frequency response of niters with a quadratic chirp for the chirp constants 
F/2ir = 0, 0.25, 0.5, 0.75, and 1, and for «L-values of (a) ir/4, (b) tt/2, (c) 3tt/4, and 

(d)TT. 
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Fig. 6 (continued). 
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which allows us to express the chirp parameter F in terms of the total 
number of periods N = L/A in the form 

F = 2tt-N-AA/A. (29) 

According to symmetry rule (b), the response of a filter with a linear 
chirp is symmetric. The discussion in Section IV implies that a reversal 
of the sign of F will not change the filter response since this operation 
simply means that we have turned the two-port structure around to 
measure the reflectivity at the other port. The filter response as ob- 
tained from a numerical evaluation of (9) using the fourth-order 
Runge-Kutta method is shown in Figs. 5a through 5d. Again, we have 
chosen for kL the values ir/4, ir/2, 3ir/4, and w. We have included the 
curves for the unchirped filters (F = 0) and shown the responses of 
five chirped filters up to a chirp constant F/2w = 3. We note that the 
chirp washes out the zeros and broadens the response. 

7.2 Quadratic chirp 

For a periodic structure with a quadratic chirp, we write the spatial 
frequency variation in the form 

AK(z) = 4>' = 12Fz*/L 3 (30) 

using F as a chirp parameter. With this we get for the phase function 

* = 4F(*/L)«. (31) 

The variation of the grating period is then given by 

AA«._a'4Y*Y, (32) 



A 2tt 



A (2z\ 
L\L)> 



and the difference between the periods in the middle of the device 
(2 = 0) and the device ends (2 = ±L/2) becomes 

AA AA(0) - AA(L/2) „ F A , Q „, 

t = a = 3 *2vr (33) 

Again, we can express the chirp parameter F in terms of the total 
number of periods by 

ZF = 2vN-AA/A. (34) 

The filter response of a periodic structure with a quadratic chirp is 
asymmetric. If we reverse the sign of F in this case, we are dealing 
with a different structure which has a different filter response. Accord- 
ing to the discussion of Section IV, this new response curve is the 
mirror image of the response of the original filter relative to the center 
frequency. The response, as obtained by a numerical evaluation of (9), 
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is shown in Figs. 6a through 6d for four values of kL. Curves are given 
for chirp parameters up to F/2t = 1. We note that with increasing 
chirps the asymmetry of the response curve increases with the side-lobe 
levels increasing on one side and decreasing on the other side of the 
center frequency. Most of the zeros in the response are preserved and 
shifted somewhat, and there is a small shift in the frequency of peak 
reflectivity. There appears to be no significant broadening of the main 
lobe of the response for the parameters selected. 

VIII. CONCLUSIONS 

We have described the filter characteristics of nonuniform periodic 
waveguides on the basis of coupled-wave theory. We have also shown 
that the filter response is described by a Riccati differential equation 
and have presented solutions for linear and quadratic tapers of the 
coupling coefficient and for linear and quadratic chirps in the grating 
period. The results can be used as an aid in filter design and to predict 
the effect of imperfections introduced during device fabrication. 
Quadratic tapers appear to be a promising choice when filter charac- 
teristics with low side-lobe levels are desired. 
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